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Part  one  of  this  report  deals  with  attempts  over  the  past  several  years  to  devise 
a  more  sophisticated  theory  of  diffusion  controlled  reactions  than  that  presented 
by  Collins  and  Kimball.  In  particular,  the  invest i gators  were  interested  in  a 
more  realistic  formulation  of  the  problem  of  high  concentration  quenching  where 
quenches  in  the  vicinity  of  the  molecule  to  be  quenched  must  be  considered.  It 
was  desired  however,  to  obtain  a  formalism  which  was  tractable  mathematically  and- 
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which  contained  parameters  which  would  be  related  to  experiment.  Part  two  deals 
with  the  photophysics  of  systems  exhibiting  molecular  association  both  in  the 
ground  and  excited  states  has  been  studied.  The  emphasis  has  been  on  kinetic 
models,  the  measurement  of  rate  constants  associated  with  these  models,  and  the 
determination  of  activation  parameters  and  equilibrium  thermodynamic  parameters 
associated  with  the  exciplex  formation  and  disappearance.  Studies  of  solvent 
effects  and  steric  effects  on  the  behaviour  of  exciplex  systems  have  been  carrie 
out.  The  case  of  rapid  equilibrium  where  the  monomer  and  exciplex  decay  with 
the  same  rate  constant  has  also  been  examined. 
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SECTION  I 


INFLUENCE  OF  DIFFUSION  ON  FLUORESCENCE 
QUENCHING. 


This  section  deals  with  our  attenpts  over  the  past  several 
years  to  devise  a  more  sophisticated  theory  of  diffusion  controlled 
reactions  than  that  presented  by  Collins  and  Kimball.  In  particular, 
we  were  interested  in  a  more  realistic  formulation  of  the  problem 
of  hif^i  concentration  quenching  where  quenches  in  the  vicinity  of 
the  molecule  to  be  quenched  must  be  considered.  It  was  desired 
however,  to  obtain  a  formalism  which  was  tractable  mathematically 
and  which  contained  parameters  which  would  be  related  to  experiment. 


1 .  S0M5  HISTORY 

Early  interest  in  diffusion  controlled  processes  centred  around 
the  coagulation  of  colloidal  particles.  These  particles  were  large 
enough  that  the  solvent  provided  an  almost  perfect  continuum  for 
diffusion.  Fick's  laws  were  invoked  and  concentration  independent 
diffusion  coefficients  employed.  It  was  assumed  that  one  could  use 
a  mutual  diffusion  coefficient  (D^  +  Dq)  to  describe  diffusion  in  a 
system  where  only  D  or  Q  was  allowed  to  move. 

Smoluchowski1  in  1917  published  a  paper  describing  the  solution 
to  Fick's  second  law  for  reactions  occurring  with  100%  efficiency  upon 
each  encounter.  Collins  and  Kimball^  modified  this  approach  to  take 
into  account  the  possibility  of  non-productive  encounters.  Their  work 
was  in  part  stimulated  by  a  paper  by  Sveshnikoff^  which  improperly 
treated  this  problem  of  inefficiency. 

Noyes^  in  a  series  of  papers  in  the  60' s  examined  the  problem  from 
the  point  of  view  of  pair  diffusion  with  essentially  the  same  results. 
Noyes'  analysis  provides  valuable  insight  into  the  origin  of  the  phen¬ 
omenon. 


2 


Others  active  in  this  fieldwere  Waite,5  Weller,6  La  Mer  and 
Unberger, 7  and  Yguerabidi,  Dillon  and  Burton.8 

9 

In  1973,  Fixman  and  Wilemski  described  a  more  general  approach 
to  the  problem.  Recently,  Andre,  Niclause  and  Ware*8  starting  with 
Fixman- Wilemski  equation  have  attempted  to  develop  more  realistic 
models  for  diffusion  controlled  and  partially  diffusion  controlled 
reactions . 

In  spite  of  all  this  theoretical  work,  very  little  experimental 
work  has  been  published  which  provides  laboratory  verification  to  the 
equations  which  are  in  common  use  for  correcting  for  transient  effects 
in  diffusion  controlled  processes.  Ware  and  Novros1*  and  Warp  and 
Nemzek*2  studied  systems  where  the  non-exponential  decay  could  be 
subjected  to  quantitative  analysis  and  a  comparison  made  between 
transient  and  steady-state  quenching.  Work  in  the  author's  laboratory 
in  collaboration  with  J.C.  Andre,  M.  Bouchy  and  M.  Niclausel3  has 
recently  resulted  in  additional  experimental  work  and  the  re-analysis 
of  some  of  the  data  from  earlier  work.  However,  the  integrated 
theoretical  output  still  exceeds  experimental  work  in  this  field 
by  a  wide  margin. 


2.  SIM>LE  THEORY 

Smoluchowski1  solved  Pick's  second  law 


If  ■  v* 


(1) 


for  the  case  where  one  particle  is  fixed  at  the  origin  and  the  other 
particle  moves  with  diffusion  coefficient  D^g.  c  is  the  concentra¬ 
tion  of  the  diffusing  species,  in  molecules  per  cm8,  and  for  radial 
diffusion 


V 


2 


(2) 


Smoluchowski  considered  every  encounter  between  particles  to  result 
in  reaction  and  therefore  used  the  boundary  conditions 


c(r,0)  -  cQ,  c(«,t)  *>  cQ,  c(o,t)  =  0 


(3) 


where  cQ  is  the  initial  concentration  of  the  diffusing  species,  r  is 
the  distance  from  the  origin,  and  o  is  the  sum  of  the  radii  of  the 
two  particles. 


The  solution  to  Eq.  (2)  with  these  boundary  conditions  is 
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c(r,t)  = 


c  1 - erfc 

o  r 


where 


erfc(x)  =  —  f  e~z  dz 

/n  J 
x 

The  flux  at  time  t,  through  the  sphere  of  radius  o  around  the  origin, 
is 


$  =  4toXd  i 

AB  ^  3r  T=a 

Taking  the  derivative  of  Eq.  (4),  one  obtains 


=  4liaDABCo  1  + 


The  observed  rate  constant  is  thus 

k(t)  *  =  4ttcjD  N '  1  +  ■  (7) 

0  i  *&) 

where  CQ  is  now  expressed  in  moles  per  liter,  and  N'  is  Avogadro's 
number  divided  by  1000. 

2 

Collins  and  Kimble  pointed  out  that  the  equations  derived  by 
Smoluchowski  contained  certain  disturbing  features.  Firstly,  at 
t  ■  0,  the  flux  as  given  by  Eq.  (6)  was  infinite.  Secondly,  the 
theory  made  no  provision  for  cases  wherein  the  probability  of  reaction 
on  encounter  was  less  than  unity.  Sveshnikoff^  had  attempted  to 
correct  for  this  second  point  by  multiplying  the  flux  by  a  factor,  a, 
less  than  unity.  Collins  and  Kinfcle  showed  that  this  procedure  was 
incorrect  since  it  implied  that  the  flux  differed  from  the  decrease 
in  the  total  amount  of  the  diffusing  species. 

Collins  and  Kinble  eliminated  these  difficulties  by  allowing  for 
the  possibility  that  not  every  particle  reacted  on  reaching  distance 
o  from  the  origin.  This  was  done  by  assuming  that  the  rate  of  reac¬ 
tion  was  proportional  to  the  probability  that  the  diffusing  particle 
was  between  0  and  (a  +  da).  That  is, 

$  *  <  4tt  0^  c(0,t)  (8) 


4' 


where  k  is  some  rate  constant.  Since  the  flux  is  also  given  by  Eq. 
(5),  one  obtains  the  condition 


c(o,t) 


AB 


9c(r,t) 

3r 


r=o 


(9) 


This  is  known  as  the  "radiation  boundary  condition",  due  to  its  use 
in  heat  conduction  problems.  Eq.  (1)  is  now  solved  with  the  boundary 
conditions 

c(r,0)  =  c 


c(®,t)  =  c 


(10) 


c(o,t)  = 


AB 


3c(r,t) 

\ 


r=o 


The  solution  with  these  boundary  conditions  is  well  known. 


c(r,t)  =  c 


,|l  -  2“~|erfc 


l/4D^t 


-  exp 


-  .  r  -  a, 

”  T“  J 


/oTIt  „ 

erfc  I 


til) 


where 


oD 


AB 


°AB  +  KO 


(12) 


Using  Eq.  (8)  and  setting  r  equal  to  a  in  Eq.  (11),  one  obtains  for 
the  flux 


4lToDABCo  r 


1  ♦ 


AB 

ICO 


1  ♦  ex  erfc(x) 
AB 


(13) 


where 
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For  x  >  1,  one  say  use  the  asymptotic  expression  for  erfc(x) 

2 

eX  *  erfc(x)  =  Cl - +  — ~r  — +  . .  0  (15) 

x/n  2x  4x  8x 

The  error  involved  in  truncating  this  expansion  is  less  than  the 
value  of  the  first  neglected  term.  For  ordinary  values  of  ic,  D^g* 
and  a,  as  will  be  shown  later,  the  term  (1/2x2)  is  equal  to  0.01  in 
approximately  10" 1®  sec.  Thus,  all  but  the  first  term  in  the  expan¬ 
sion  can  be  neglected,  and 


x  1 

•  erfc(x)  =  - 

x/n 

The  flux  then  becomes 


(16) 


<&  =  4irD.0a  c 
AB  o 


1  + 


(17) 


where 


1  + 


AB 

KCT 


(18) 


The  time  dependent  rate  constant,  k(t),  is  then  given  by 


k(t)  =  ±  =  4irD^gdN' 
o 


1  + 


(19) 


This  equation  is  of  the  same  form  as  Smoluchowski's  equation,  Eq. 

(7),  except  that  o  has  been  replaced  by  a. 

It  is  of  interest  to  examine  the  flux  as  the  time  approaches 
zero.  For  small  values  of  t,  one  cannot  use  the  asymptotic  expansion 
of  erfc(x),  but  instead  must  integrate  the  Taylor's  series  expansion 
of  e-*2  to  obtain 

2  f  _r2  2x  x2  x^ 

erfc(x)  *  —  J  e  Ax  =  1  -  ^  (1  -  ~  *  *  0  (20) 

x 

2 

x 

On  multiplying  this  expression  by  the  Taylor's  series  expansion  of  e  , 
one  finds  that  the  flux  at  small  t  is  given  by 
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4iraDABCo 


D 


1  + 


AB 

<a 


1  + 


K0 

dab 


Cl 


(21) 


The  flux  as  given  by  this  equation  is  finite  as  t  approaches  0.  At 
t  =  0,  the  flux  is  equal  to  4itct2kco. 

Eq.  (19)  is  the  basic  equation  which  is  used  in  the  analysis  of 
the  fluorescence  quenching  data.  The  derivation  of  this  equation 
depends  critically  on  the  validity  of  truncating  the  asynptotic  expan¬ 
sion  of  ex^  erfc(x) .  This  truncation  neglects  the  term  (2X2)-!  where 


x 


(°AB  +  KO) 


CD 


AB 


(22) 


2-1 

In  order  to  truncate  the  expansion  in  this  fashion,  (2x  )  must  be 
much  less  than  unity.  It  is  desirable  to  determine  the  range  of  a, 
D^jj,  and  t  over  which  Eq.  (19)  can  be  used.  To  do  this,  one  must 
estimate  the  value  of  k,  and  then  calculate  values  of  (2x2)"!  for 
given  a,  D^,  and  t. 

As  can  be  seen  from  Eq.  (8),  k  has  the  units  of  cm  sec-1. 
Noyes4  assumed  that  one  can  relate  K  to  the  rate  constant,  k„,  cal¬ 
culated  from  the  kinetic  theory  of  gases,  using  the  relation* 


2 

K  =  4iro  tc 

g 


(23) 


where  Kg  is  on  the  order  of  10 


-10 


cm^  molecule  1 


-1 


sec  .  Let  us  take 

5  x  10~8  cm  as  a  reasonable  estimate  of  0.  Then  in  non- viscous 
media,  with  D^b  =  5  x  10-5  cm2  sec*l,  x  is  equal  to  5.9  x  105  v't. 
When  t  is  10  psec,  (2x2)-^  is  found  to  be  0.14,  which  is  not  neglig¬ 
ible  compared  to  unity.  When  t  is  100  psec,  (2x2)-1  is  equal  to 
0.014,  which  is  small  compared  to  unity. 


“6  2  “1 

In  more  viscous  media,  with  D  =  5  x  10"  cm  sec  ,  x  is  equal 
to  1.47  x  106  /t.  When  t  is  10  psec,  (2x^)”l  is  0.023.  When  t  is 
100  psec,  (2x2)-*  is  0.0023. 


It  is  thus  apparent,  that  even  if  the  estimate  of  k  is  only 
approximately  correct,  the  truncation  of  the  expansion  is  valid  for 
nsec  time  scales  in  viscous  media. 


3.  APPLICATION  TO  FLUORESCENCE  QUENCHING 

Consider  the  following  fluorescence  reaction  mechanism: 


7 


A  A* 
kl 

A*  — A 
k2 

A*  — =-*■  A  +  hv 

A*  +  B  ■  k ^ — *■  Products 

The  rate  constant  k(t)  is  given  by  Eq.  (19),  and  (kj  +  k2)_1  is  equal 
to  t0,  the  unquenched  lifetime  of  the  system.  For  6-pulse  excitation 


ilill  .  [A*]  (-  ±  -  k(t)[B])  (24) 

o 

When  the  concentration  of  B  is  much  larger  than  [A*]0,  then  B  can  be 
considered  to  be  time  independent,  and 

tpLr  =  exp  (-at  -  2b /t  )  (25) 

Lo  J 

where  [Aq*]  is  the  concentration  of  A*  at  time  0,  and 

a  -  •—  +  4ird  D  N’  [B]  (26) 

o  ^ 

b  =  4/irD^  a2  N*  [b]  (27) 

The  decay  law,  G(t) ,  is  predicted  to  be 

G(t)  =  exp  (-at  -  2b/t  )  (28) 


The  response  of  the  system  under  steady  state  illumination  is 
obtained  by  considering  excitation  by  a  step  function  F0.  If  the 
number  of  excited  molecules  of  A  is  given  by  nA(t) ,  we  have  that 

t  t 

nA(t)  =  I  F(t  -  t')  G(t')  dt’  =  Fq  J  G(t')  dt' 

0  0 

when  F(t)  is  a  step  fimction.  On  obtaining  an  expression  for  n^(t), 
one  then  lets  t  -*•  00  in  order  to  arrive  at  the  number  of  molecules, 
"A(“)»  in  the  steady  state. 
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I 


\ 


The  result  is 


12 


I 

o 

I 


nA^B=0 

nAM 


(1  ♦  4todaqn'£b]t0) 


-a/! 


exp 


erfc 


- 

s* 


I 

(29) 


This  expression  may  be  used  for  analysis  of  steady  state  data  obtained 
under  conditions  where  the  time  dependence  of  the  rate  constant  is 
expected  to  be  important. 

Eqs.  (28)  and  (29)  are  frequently  used  for  the  analysis  of  dif¬ 
fusion  controlled  reactions. 

Experimental  studies  of  the  validity  of  Eqs.  (28)  and  (29)  indi¬ 
cate  the  following: 11 » 12 , 14, 15 

(a)  o  is  always  greater  than  the  collisional  distance  o. 

(b)  Calculated  values  of  a  from  steady  state  measurements  are 
greater  than  values  of  o  from  decay  law  analysis;  with  dramatic  depart¬ 
ures  in  some  cases,  of  the  Stem-Volmer  plot  from  that  predicted  from 
the  decay  law  parameters. 


4.  AN  IMPROVED  MODEL 

A  more  realistic,  but  still  simple  model,  has  recently  been 
proposed  by  Andre,  Niclaus  and  Ware.*0 

Theoretically,  it  is  possible  to  determine  variations  of 
apparent  rate  constant  of  reaction  ka(t)  versus  time  by  resolution  of 
the  classical  equation  of  Wilemski  and  Fixman:9 

3<}>(r,t)/3t  +£<)>(r,t)  =  -k(r)<f>(r,t) ,  (30) 

where  <j>(r, t)  is  the  ratio  of  the  configuration  distribution  function 
of  the  molecules  B  to  the  mean  concentration  of  B,  is  a  time  evolu¬ 
tion  operator  relating  to  the  diffusion  processl6 

£x  =  -DV2x  +  KVY  •  Vx  +  KV^x,  (31) 

with:  V  the  electrical  potential  around  the  molecule  A*,  K  a  mobility 
coefficient. 

k(r)  is  a  rate  constant  (expressed  in  s  *)  of  reaction  of  an  A* 
molecule  with  a  B  molecule  at  the  distance  r. 

The  analytical  or  numerical  variations  of  k(r)  versus  r  have  not 


9 


been  reported  in  the  literature  (except  for  Forster  energy  trans¬ 
fer)  .  Except  in  the  case  of  chemical  reaction  between  ions,  the 
value  of  ¥  is  not  known.  However,  one  can  assume  that,  in  the  case 
of  uncharged  molecules,  ¥  takes  a  significant  value  only  for  values 
of  r  near  of  o. 18 

The  model  of  Andre,  Niclause  and  Ware  is  a  diffusion  model  con¬ 
sisting  of  two  distinct  volumes  in  the  reacting  system  centered  on  a 
molecule  A*  towards  which  molecules  B  diffuse: 

(i)  a  reaction  range  between  o  and  o'  wherein  the  chemical 
reaction  occurs  with  a  chemical  rate  constant  practically  infinite; 
we  assume  that  no  diffusion  occurs  in  this  range,  i.e.,  <Kr,t)  =  0 
for  r  £  [a, a* ], 

(ii)  an  external  volume  (r  >  a')  in  which  no  chemical  reaction 
but  only  diffusion  occurs. 

Then,  we  have  to  solve  the  following  system: 

for  r  €  [a, a' ]  and  t  >  0 

<Kr,t)  =  0; 

for  r  >  a' 

9<|>(r,t)/3t  +J:<H*,t)  =  0, 

with 

J:<K  r,t)  =  -DV2<p(r,t).  (32) 

Under  these  conditions,  the  value  of  the  apparent  rate  constant 

1  IQ 

is  given: 1,13 

ka(t)  =  4TrNa’D[9a(r,  t)/9r]0, , 
which  becomes 

ka(t)  =  kQ[l  +  0  '  (TTDt)  ”^] ,  (33) 

where  k  =  4tiN0'D. 
o 

If  A*  is  a  singlet  excited  state,  one  can  measure  the  intensity 
of  fluorescence  of  A*  in  the  presence  (I)  or  in  the  absence 
(I0)  of  B,  all  other  experimental  conditions  being  constant. 

Theoretically,  one  determines  the  variations  of  I  with  concentra¬ 
tion  [b]  taking  into  account  that  for  t  =  0,  there  is  a  probability  P 
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of  finding  a  molecule  B  inside  the  chemical  reaction  volume  Vj  [Vj  = 
4/3  it  (o' $  -  a3)].  Then  one  can  postulate 

(i)  practically  instantaneous  deactivation  of  A*  having  at 
least  one  B  molecule  in  V^, 

(ii)  deactivation  from  B  molecules  reaching  a'  after  diffusion 
(Smoluchowski) . 

All  other  conditions  being  constant,  for  the  steady-state 
experiment,  the  ratio  of  I  over  I  can  be  expressed  by: 

00  /  oo 

T  ■/<*•>[: B]Wj  (34> 

o  /  0 

On  the  contrary,  by  flash  excitation  where  one  measures  the 
relative  variations  of  A*  versus  time,  only  dynamic  quenching  will  be 
observed. 

20 

The  generalized  distribution  of  Smoluchowski  gives  the  proba¬ 
bility  of  having  a  given  distribution  of  B  molecules  between  o  and  o' 
around  A*. 

If  [B]  is  the  average  concentration  of  B,  N  the  Avogadro's 
number,  the  probability  P(n)  to  have  n  molecules  B  in  Vj  is: 

-NV.  [B]  _ 

P(n)  =  e  CNVj[B])  /n!,  (35) 


or 


P(n)  =  e  Wwn/n! 


with 


w  =  NVj[B]. 

The  probability  P  to  have  at  least  one  B  in  is  expressed  by: 

oo 

_  „  -w  n ,  ,  ,  -w 

P=E  e  w/n!=l-e 

i=l 

Thenj.  if  (A)0  molecules  are  excited  at  t  =  0,  we  obtain 
(A)0e_NVlLBJ  excited  molecules  after  this  very  fast  deactivation 
process. 


The  following  scheme  for  the  disappearance  of  A*  can  be  used 


A*  - ^ — *■  . . .  unimolecular  relaxation  to  ground  state 

kaoo 

A*  +  B  - ►  ...  bimolecular  quenching 

then 

d(A*)/dt  =  -  [k  +  ka(t)[B]](A*) 
leading  to 

(A*)(t)  =  (A*)Qexp[-NV1[B]]  x  exp{-  [k  +  kQ[B] 

+  20I(TrDt)'34ko[B]]t},  (36) 

and 


00 

'iM  / 


exp[NV.  [B]]  I  exp(-kt)dt 


I  ® 


/ 


-it. 


exp{-  [k  +  kQ[B]  +  2a' (irDt)”  ■2kQ[B]]t}dt 


then 


!  NVjCb] 

f  '  TOJ-  [l  *  k°t[B]]' 


where 


X  = 


2  V.W  o' 


T7  TT  > 


IT*  [1  ♦  k0T0[B]T*  (Dlo) 


and 


A 


1(A)  =  J  exp(-u  -  Au*)du 
0 


(37) 


(38) 
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■  1  -  i  xAxp(j  X2)Erfc(y  X). 

Then,  by  flash  excitation,  one  measures  the  relative  variations 
of  (A*)(t)  (Eq.  (36))  due  to  unimolecular  relaxation  and  to  the  reac¬ 
tion  of  A*  and  B  molecules  reaching  a'  (the  true  transient  effect 
occurs  in  o'  (iTDt)"^) .  On  the  contrary,  by  steady-state  measurement, 
one  other  effect,  due  to  the  reaction  of  A*  with  B  presents  in  is 
su^e|,i^posed  on  the  other  (apparent  transient  effect  due  to 

e  1  J) . 

For  small  quencher  concentrations,  the  first  and  the  second 
terms  of  the  Taylor  series  of  IQ/I  versus  [B]  are: 

I  , 

=  1  +  {NVj  +  kQTo[l  +  a'(Dxo)'^]}[B]  +  ...  (39) 

All  other  terms  in  this  expansion  can  be  neglected.  Then,  the 
apparent  constant  of  quenching  Kgy  is  equal  to: 

KSV  =  NV1  +  koToC(1  +  (40) 

„  For  most  of  the  studies  which  are  reported,  values  of  o'  around 
10  A  (j*  1)  were  obtained  (for  example,  refs.  11,12,14,15).  But,  the 
calculated  collision  value  is  always  around  6  X.  Then,  under  these 
conditions,  the  importance  of  is  about  NVj  *  1  mole-1  £. 

For  low  viscosity  solutions  and  for  excited  molecules  with  life¬ 
times  greater  than  10  ns,  one  can  neglect  this  static  effect.  On  the 
contrary,  in  order  to  observe  transient  effects,  one  needs  a  value  of 
a* (Dt0)-*5  as  great  as  possible  which  can  be  obtained  by  lowering  D 
(or  increasing  viscosity  n)  and  by  decreasing  k0,  i.e.,  by  increasing 
relative  importance  of  static  quenching  term  (see  Eqs.  (38)  and  (39)). 


5.  A  HI  (HER  LEVEL  OF  SOPHISTICATION 

A  more  rigorous  theory,  but  one  that  still  yields  relatively 
simple  equations,  has  been  developed  by  Andre,  Bouchy,  Niclause  and 
Ware.13 

One  again  starts  with 

♦^(r.t)  =  -k(r)  •  <Kr,t)  (30) 

Knowing  £and  k(r),  it  is  possible  to  determine  the  variation  of 
$(r,t)  with  t,  by  using  the  Laplace  transform  or  numerical  tech- 


13 


niques.1^  The  apparent  rate  constant  is  given  in  M_**sec~*  by  the 
expression 

o' 

ka(t)  =  J  4nr2Nk(r)<j>(r,t)  dr 
a 


In  the  case  where  a  chemical  reaction  occurs,  the  variation  of 
k(r)  is  not  well  known  and  a  simplified  model  must  still  be  enployed. 
Reaction  or  interaction  between  molecules  requires  the  overlap  of  the 
external  molecular  orbitals  of  the  two  molecules.  The  value  of  the 
overlap  of  molecular  orbitals  is  known:  it  is  a  decreasing  function 
of  r  and,  beyond  a  value  of  r  of  some  A,  there  is  practically  no 
overlap. 

In  our  model,  we  take  k(r)  to  be  practically  constant  and  equal 
to  some  value  K  for  distances  smaller  than  o',  i.e.,  for  a  <  r  <  o' . 
At  distances  greater  than  o',  k(r)  is  set  equal  to  zero.  Thus,  the 
equation  for  K(t)  becomes 

a' 

Ka(t)  =  47TNk  J  $(r,t)r2dr  (41) 

0 

We  have  proposed  a  diffusion  model  considering  of  two  distinct 
volumes  in  the  reacting  system  centered  on  a  molecule  A  toward  which 
molecules  B  diffuse: 


-  a  reaction  zone  consisting  of  the  volume  between  a  sphere  of 
radius  o'  and  a  sphere  of  radius  0  wherein  the  chemical  reaction 
occurs  with  a  rate  constant  k;  we  assume  that  no  diffusion  occurs  in 
this  range,  i.e.,  <f>(r,t)  is  independent  of  r. 

-  an  external  volume  (r  >  o')  in  which  no  chemical  reaction  but 
only  diffusion  occurs. 

The  value  of  (?'  and  the  rate  constant  k  are  determined  by  the 
variation  of  ip  and  by  the  particular  value  k(r)  in  the  vicinity  to 
the  molecule  A.  The  model  is  illustrated  in  Fig.  1. 


For  r  >_  a’,  taking  ip  =  0  for  uncharged  particles. 


||  =  D 72<p  =  D 


+  i 

l  3r2  r  3r 


(42) 


For  o  £  r  <  o’,  the  flow  of  molecules  B  at  the  distance  o'  of  A 
has  to  be  equal  to  both  the  accunulation  of  B  inside  the  reaction 
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volume  and  the  disappearance  of  B  in  the  chemical  reaction  A  +  B  ■+■  C. 
Thus,  we  must  have 


4ttN0'2D| 


3r 


r,<[B]>=iuN(a'3-o3)(||Ja(<[B]> 

♦  jttNCct'3  -  a3)  •  kC<j>)0,<[B]> 


(43) 


The  macroscopic  bimolecular  rate  constant  kc  (in  mole”  s”  )  is 
related  to  the  rate  constant  k  by  the  equation 


o' 

^  =  kc[A][B]Gf  =  [A][B]a,N  J  4irr2kdr 

a 

then 

kc  -  ~  ttN(o' 3  -  o3)k  «  NVjk  (44) 

By  including  this  macroscopic  term  in  Eq.  (43),  one  obtains  the 
following  relation: 


koa‘(lf)0-  ■  Nvi(f)c  *  kc<«0- 

where  k  *=  4nNa'D. 
o 

Eqs.  42  and  45  together  define  the  systems. 


(45) 


It  is  convenient  to  use  the  dimensionless  terms  p  *  r/o'  and 
t  ■  Dt/o'2;  eqs.  42  and  45  become 


f 

1 

f  0  1 

■ 

3 

9di 

(ip) 

1  ”3 

1  *  j 

m 

■ 

[A 

(46) 


techniques 


This  system  lends  itself  to  resolution  by  Laplace 


The  transform  y(p)  of  $(t) ,  is  defined  as 


rife 
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00 

y(p)  3  J  e'pt4>(t)dt 


Then,  we  have  (cf.  Ref.  13) 


[y(p)]i  =  yi(p)  B  4  1  ~  - ~ - 

FL  ap  +  *^p  +  B  +  lJ 


1  -  in 


6  =  kc/ko 


The  value  4>i(t)  of  the  distribution  function  at  the  distance 
o'(p  3  1)  gives  several  analytical  solutions  which  depend  on  the 
value  of  the  quantity  a(8  ♦  1) . 


(i)  Case  Where  q(g  +  1)  =  0 

♦l(T)  3  +  g-~  eXp  C(3  +  1)2t]  Erfc  [CP  ♦  1)*^  3  (48) 

11  12 

This  relation  is  similar  to  that  proposed  by  Ware  et  al.,  ’  by 
Owen21  and  Collins  and  Kinfcall.2  However,  this  equation  corresponds 
to  the  case  where  o'  =  a,  whereas  these  authors  assume  a'  >  a,  but  do 
not  take  into  account  the  accumulation  of  B  corresponding  to  the  term 

NV^at/at) . 


(ii)  Case  Where  a(B  +  1)€]0,  1/4 1 


a  „  1  2a8 

6  4  1  ~  Y(1  + 


+  ref?  y)  exp 


with  y  3  /I  -  4(6  +  l)a  . 


ex>,[(i*H  t1  Erfc 

’,[(t5J1)  ,]Erfc[(i2jJ:)/r] 
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(iii)  Case  Where  a(8  +  1)  a  1/4 


^(x)  -  (1  -  4a2)  -  4a2 


2a 


exp  — j  Erfc 
4a 


2a, 


-  —  A/tt 
a  ' 


(50) 


(iv)  Case  Where  g(g  +  1)  >  1/4 


1  .  2ia8 

i  ! 

[fi  -  tr! 

2' 

r  1  -  iy' 

e  +  i  y 

1  +  iy'  exP 

Ll  2a  i 

I  T 

crxci  | 

U 

l  2a 
i 

exp 


1  +  iy1 


U 


2a 


2  1 
T 


1  -  iy' 

with  i2  =  -1,  and  y 1  =  /4a (8  +'  1)  -  1  . 
From  above,  one  obtains 
ka(r)  *  NV^fx)  =  k^fx) 


Erfc 


1  +  iy' 
2a 


/r 


/T 


(51) 


(52) 


The  variations  of  k  (x)  are  therefore  given  by  the  expression  of 

^(T). 

At  large  values  of  time,  the  above  relations  can  be  simplified 
by  using  the  Taylor  series  expansion  of  the  transform  yj(p).  One 
obtains 


1 


3  rr i 


1 


(53) 


and 


k  k 

*  ittV 


1  + 


k  +  k  t 
c  o  /irDt 


(54) 


This  simplified  relation  is  identical  to  that  proposed  by  Nemzek 
and  Ware. ^ 

The  limit  of  k  (t)  as  t  approaches  infinity  is 

a 


Lim  k  ■  •: - ; — rr— 

t-  *  1  *  kc/ko 


or 


. .  1  1  1 
Lim  F“  =  T  +  IT 

t-*»  a  c  Ko 


and  when  kc/kQ  approaches  infinity,  k&(t)  becomes 


-  k0  1  + 

/irDt 


which  is  identical  toEq.  (19),  except  for  the  value  of  0  which  has 
been  replaced  by  o'. 

In  the  analysis  of  experiments  on  fluorescence  quenching  using 
conventional  flash  excitation,  Eq.  (54)  can  be  used.*^  if  the  pheno¬ 
menon  is  recorded  during  a  very  short  interval  (several  picoseconds) , 
the  former  expression  of  ka(t)  is  no  longer  valid.  A  Taylor  series 
of  Yi(p)  in  the  term  l/t^>  leads  to  the  following  expression  of  k  (t) 
valid  for  short  times  a 

• 

k.  (t)  ■  k  « w. »»  -  t2 


3/fr  a* 


15a4/rr 


t5/2  + 


Time  evolution  of  the  apparent  rate  constant  ka  is  given  in  Fig. 
2  by  using  the  three  equations  51,  54  (for  long  times)  and  56  (for 
short  times).  The  range  of  validity  for  each  equation  is  clear  from 
the  diagram. 


6.  CONSEQUENCES 

Consider  the  following  scheme  for  the  deactivation  of  A*: 

A,. 

.  -  A  ♦  h\>„ 


A  or  T 


A*  +  B 


Vjft) 


After  6-pulse  excitation, 

■  "  t(k  ♦  1jc[b]01  (t)  ][A*] 
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Transforming  to  non-dimensional  variables  as  before, 

,.2 

c,  «v  w 
C 

_ _ —  *  - 

D 


, 1 ,  , 
Ta*J  d 


.  2  k  a 

kl3  .-V^w 


Then,  if  we  define 


Z(t)  =  log 


tat; 


*0 


(58) 


and 


ka 


,2 


ko'2 

a2  =  -V  w 


we  obtain 


dZ 


£  -  -Ai  -  Vl<T> 


(59) 


It  is  possible  to  resolve  this  system  using  Laplace  transforms.  One 
obtains 

^1 

PZ(p)  -  Z(0)  =  -p  -  A^jCp) 

with 


yl(p)  sf  e_PtMt)dt 

0 


leading  in  Laplace  transform  space  to 


Z(P) 


_ 1 _ 

ap  +  /p  *  &  *  1 


(60) 


It  is  also  possible  to  obtain  the  value  Z(x),  of  the  decay  func¬ 
tion  of  A*.  Several  analytical  solutions  are  obtained  which  depend 
on  the  value  of  the  product  a(&  ♦  1)  (cf.  Ref.  13): 
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(i)  Case  Where  a(8  +  1)  c  0 

A„ 


Z(t) 


A,  + 


1  6  +  1 


2A-8 

T - 1 U-/F7FF 

ce  + 1  r 


A-8  ~ 

+  - — =•  [1  -  expect  +  1)  T)Erfc(C$  +  1)/F  )] 

(6  +  1 V 


(61) 


(ii)  Case  Where  a(8  +  1)€]0,  1/4 C 

2A_g 


Z(T)  -  - 


A,  + 


a2  \ 


1  3+1 


O  M  _ _  '''O' 

2  /FA  2 


A,8(l  -  a(6  +  D) 


8tt3A2B 


(1  +  Y) 


(3  +  1) 


1  4a2r 

- 3  exp  - 


(6  +  1) 


(1  -  Y) 
with  y  =  A  -  4a(&  +  1) 


1  (  4a2x 

- =  exp  - : 

1(1-  y) 


(1  +  Y) 
Erfc1 


Erfc 


2a 


1  +  Y 


ft 


2a 


1  -  Y 


ft 


(62) 


(iii)  Case  Where  a(3  +  1)  =  1/4 

A„ 


Z(T)  =  - 


A,  + 


1  8  +  1 


a28 


2A  6 

"  ftft 


3A28 


(3  +  1)' 


4(8  +  ir  L 


2(0  +  1)  2x  -  ||  Erfc 


4(8  +  1)' 
I  rr\ 


ft 

l2ai 


•  exp] 


4a 


(63) 


(iv)  Case  Where  a(0  +  1)  >  1/4 
I  A_ 


Z(T)  *-  +^11- 


8a3A  8 
♦  - : -  1 


2A-B  A,8(l  -  a(0  +  1)) 

- r  /t/  tt  ♦  — - ; - 


(6  +  1)' 


(8  +  1)' 


Y' 


— 

i 

[  4a2  T1 

Erfc 

2a  ft 

t  exp 

(1  -  iY') 

L  (1  -  iY')2  J 

L  (1  -  iY’)  J 

MM 


4 


(1  +  iy')' 


exp 


4a 


L  ci  +  iy')2 


Erfc 


2a 


ci  +  iy') 


/r 


(64) 


with  =  -1  and  y'  =  /4a(0  +  1)  -  1 

These  different  expressions  lead  to  variations  of  the  value 
A„ 


y  s  z(x)  + 


A1  +rrr 


X  + 


2A -g  __ 
/t/tt 


(3  +  1)' 


as  it  is  shown  in  Fig.  3. 


As  can  be  seen  from  Fig.  3,  for  values  of  t  of  the  order  of  mag¬ 
nitude  of  unity,  or  greater  than  unity,  the  difference  between  Z(t) 
and  the  decay  function  Z'(t),  which  contains  only  terms  proportional 
to  t  and  /F,  i .  e . , 


Z'(t) 


+ 


A2  1 
6  +  ljT 


2A,0 

- _  /?7? 

(3  ♦  1) 


is  practically  constant. 

The  expression  of  Z'(t)  can  be  obtained  (Eq.  (54))  from  the  simp¬ 
lified  expression  of  4>l(T).  then  1  (x)  .  The  curves  presented  in 
Fig.  3  show  only  the  difference  between  Z  and  Z'  divided  by  A2> 

(i)  Z(x)  £  Z'(x)  for  a($  +  1)  ^  1  when  <J>i'(x)  is  always  greater 
than  <^i(t).  This  occurs  particularly  in  the  case  presented  in  paper 
I;1  the  expression  for  Z(x)  in  this  important  case,  when  the  reaction 
is  practically  limited  by  diffusion,  is  given  by 

Z(t)  =  -NV^B]  -  [(k  +  k0[B])t  +  2koBa«  /(t/irD)  ] 


(ii)  Z(T)  >  Z'(t)  for  a(0  +  1)  <  1.  This  result  is  obtained 
for  the  conditions  where  <J>i'(x)  can  be  less  than  (x) .  The  results 
can  be  understood  with  the  aid  of  the  qualitative  curves  in  Fig.  4 
and  from  the  Taylor  series  of  Z(x)  obtained  for  large  values  of  x. 


Z(x)  =  - 


A1  +  3  +  1 


2A20 


A„3 


A20[2a(0  +  1) 


(3+1)' 

«  1 


/t/t\ 


C3  +  1)' 


[a(0  ♦  1)  -  1] 


(3  ♦  1) 


/frT 


\ 


(65) 
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As  can  be  seen  in  Fig.  3,  for  small  values  of  r  (^0.2),  termin¬ 
ation  of  the  expression  (65)  after  the  constant  term  yields  a  good 
approximation  of  the  time  evolution  of  Zfr).  By  replacing  A^  and  A., 
we  obtain 


Z(t)  * 


k  k 

k  *  F^nr  « 

C  0 


k  \ 

t - [B]  -  NV'[B]  (66) 

(k  +  k  y  /irDt 


i.e.. 


Z(t)  =  -at  -  b/F  - 


with 


NVj' 


=  4irNa 


k„2k 

,3  co 


o.  *  o3  LI  loi 


C  0 


afl  ♦  -  1 


Using  the  approximate  form  of  the  above  equations,  it  is  then 
possible  to  analyze  data  from  the  following  experiments: 

(i)  Single  Photon  Counting.  One  can  obtain  the  values 

k 


k  k 
o  c 

k  +  k 
c  o 


and 


k  +  k 


o  S 


Also,  the  variation  of  kc/k0  can  be  studied  by  measuring  the  effect  of 
viscosity  to  the  reaction  between  A*  and  B.  Then,  knowing  kD,  kc  can 
be  calculated.  It  is  thus  possible  to  calculate  all  the  essential 
kinetic  parameters  of  the  reaction. 

(ii)  Steady-State  Fluorescence  Measurements.  When  the  lifetime 
of  A*,  in  the  absence  of  B,  is  short  (less  than  10  ns)  it  is  not  pos¬ 
sible  to  find  a  good  decay  function  of  the  form  exp(-at  -  b/F  -  c)  in 
which  a,  b  and  c  are  constants.  For  these  conditions,  only  steady- 
state  fluorescence  measurements  lead  to  the  determination  of  the  kin¬ 
etics  parameters  for  the  reactions  between  A*  and  B. 

When  A*  is  a  singlet  excited  state,  we  can  measure  the  quantum 
yield  of  fluorescence  of  A*  in  the  presence  (I)  or  in  the  absence 
(Iq)  of  B,  all  other  experimental  conditions  being  constant. 

Then,  all  other  conditions  being  constant,  for  the  steady-state 
experiment,  the  ratio  of  i  over  I  can  be  expressed  as  done  above  by 
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with  tq  =  1/k ;  with  I  (A')  =  J  exp  (-a  -  A'/u  )du  (1)  and 


.  ,  k  2k  t  [B] 
-  2q  c  o  o 

(kc  +  kQ)2 


/  k  k 

("jr? rToW 

c  o 


For  small  quencher  concentrations,  the  first  and  the  second 


terms  of  the  Taylor  series  of  versus  B  are: 


T  1  + 


-r-  «  1  +  <  NV  '  ♦ 


1  k  +  k  o  k  +  k 


o  c 


0  c  /Dr 


[B]  ♦  .. 


All  other  terms  in  this  expression  can  be  neglected.  Then,  the 
apparent  constant  of  quenching  Kgy  is  equal  to: 


Ksv =  "V  ♦  rrrTo 1  *  r~h 

o  c  o 


o +  Kc  sr 

o 


For  short  times  (smaller  than  about  1  ns),  a  Taylor  series  of 
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Z(p)  lead  to  the  following  expression  of  Z(t): 


Z(T) 


(Aj  +  A2)T  + 


A26  t2 


A2g  '  8  t5/2 
a2  15/if 


A20  t3 

+  -^-Cl  -a(B  +  1))  ... 


(69) 


7.  EXPERIMENTAL  RESULTS 

In  their  study  of  the  quenching  of  1,2-benzanthracene  by  CBr4, 
Ware  and  Nemzekl2  verified  that  the  decay  law  for  this  system  when 
propane  diol  was  used  as  solvent  was  not  a  single  exponential,  but 
could  be  adequately  fit  by  exp{-at  -  bt3*}.  This  is  illustrated  in 
Fig.  4.  From  Eqs.  (26)  and  (27),  the  a  values  were  obtained  along 
with  the  values  of  D.  These  are  listed  in  Table  1.  If  these  values 
are  used  in  Eq.  (29),  they  fail  to  reproduce  the  steady  state  beha¬ 
viour  by  a  wide  margin  (see  dashed  line.  Fig.  5).  On  the  other  hand, 
the  diffusion  coefficients  are  reasonable. 


This  discrepancy  is  resolved  by  using  Eq.  (38)  and  seeking  the 
best  values  of  o'  and  D  consistent  with  the  set  of  data  represented 
by  the  points  in  Fig.  5.  a  was  determined  by  Lebas '  theory. 22  The 
values  of  o'  and  D  are  in  excellent  agreement  with  those  obtained  from 
exp{-at  -  bt1*}.  This  is  illustratedoby  the  solid  lines  in  Fig.  5 
which  were  calculated  from  o'  =  8.8  A  and  the  values  of  D. 


The  decay  law  exp{-at  -  bt^}  has  also  been  verified  for  the 
quenching  of  pyrene  by  biacetyl. 23  From  Eq.  (66),  the  values  of  a 
and  b  yield  [kc/(k0  +  kc)]a'  and  D.  This  approach,  when  used  on 
decay  data  analyzed  either  by  iterative  reconvolution  or  deconvolution 
via  the  fast  Fourier  transform  technique, 24  yielded  values  of 
[kc/(kc  +  k0)]o'  of  9.1  A  and  D  =  1.53  x  10-6  cm2/sec  for  the  quench¬ 
ing  of  pyrene  by  biacetyl  in  75%  cyclohexane  -  25%  cyclohexane  sol¬ 
vent  at  22°  C  (viscosity  =  42  cp) .  The  fit  is  illustrated  in  Figs. 

6  and  7.  Eq.  (66)  requires  the  terms  a  and  b  to  depend  on  concentra¬ 
tion.  This  is  illustrated  in  Fig.  8.  From  the  variation  of  the 
quenching  with  viscosity,  kc  is  estimated  to  be  about  10*0  M_1  sec-*. 
At  high  vigcosity  where  the  value  of  [kc/(kc  +  k0)]o'  was  determined 
to  be  9.1  A,  kc  >>  kD  and  thus  a 'o  -  9.1  A.  From  the  molecular  dimen¬ 
sions,  22  one  calculates  0  =  6.5  A  in  keeping  with  the  general  pattern 
observed  in  diffusion  controlled  processes. 


A  further  example,  where  the  formalism  presented  above  has  been 
employed,  is  the  system  of  POPOP  quenched  by  CBr4  in  mixed  cyclohex- 
anol- cyclohexane  solvents. 23  As  can  be  seen  from  Fig.  9,  strongly 
non-linear  Stem-Volmer  plots  are  obtained.  Analysis  using  Eq.  (68) 
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yielded  the  [kc/(kc  +  ko)]0'  =  9  A  over  a  range  of  viscosity  of  1.3 
to  24  cp.  From  the  analysis  of  the  variation  of  [(kQ  +  kc)/kc]o'  as 
a  function  of  diffusion  coefficient,  kc  is  estimated  to  be  3  x  10*° 

W1  sec~l.  For  the  viscosity  range  of  the  experiment  of  Ware  and 

Nemzek*2 

0.98  <  k  /(k  +  k  )  <  0.998 

—  c  v  c  o 

Thus,  POPOP  quenched  by  CBr4  in  propane  diol  should  be  essentially 
diffusion  controlled. 

12 

This  lends  support  to  the  assumption  of  Ware  and  Nemzek  that 
the  quenching  of  1, 2-benzanthrancene  by  CBr4  is 
controlled  in  propane  diol. 

Application  of  the  above  formalism  to  fluorescence  polarization 
studies  of  quenching  have  also  been  made  by  Andre,  Bouchy  and  Ware. 23 
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TABLE  1.  Values  of  c*  and  D  which  give  the  best  fits  to  decay  curves  1 

(from  Nemzek  and  Ware)12  and  to  steady-state  experiment  of  I 

1,2-benzanthracene  quenched  by  CBr^  in  propane-diol . 


o'  (A)  D  106  cm2  s-1) 

« °r  ic°r  r 


CBr^  Concentration  (M) 

Tst 

25  UC 

35*0 

15  °C 

25  DC 

35°C 

from  decay  curves 

0.098 

9.0 

8.5 

8.5 

0.23 

0.50 

0.88 

0.18 

9.0 

9.0 

8.0 

0.24 

0.44 

0.95 

0.29 

8.0 

7.0 

7.0 

0.28 

0.66 

1.1 

steady-state  experiments 

(classical  model) 

(Ref. 

12) 

11.2 

- 

- 

0.287 

- 

- 

steady-state  experiments 

(Ref. 

23) 

8.8 

8.8 

8.8 

0.37 

0.59 

1.01 

* 


Figure  1.  Separation  of  space  into  two  volumes  where  chemical  reac¬ 
tion  (r  <  o')  and  diffusion  (r  >  o')  occur.  (A)  Distribution  at 
t  *»  0.  (B)  Distribution  at  t  >  0. 

Figure  2.  Substitution  of  the  true  value  of  4>i(t)  by  Taylor  series 
(a  =  0.3  and  3  =  0.5):  (1)  True  analytical  function  of  <|>i(x); 

(2)  Taylor  series  for  short  values  ofx(Eq.  56);  (3)  Taylor  series 
for  large  values  of  t  (Eq.  54). 

Figure  3.  Variation  of  Y  =  {Z(x)  +  [Aj  +  A2/(B  +  1) ]x  +  2A2B/(S  ♦ 

1)2  (x/ir)5S}/A2  versus  x  (see  text).  3=5.  (1)  a  =  0;  (2)  a  =  1/6; 

(3)  a  =  2/3;  (4)  a  =  1/3. 

Figure  4.  Open  points:  flashlanp.  Solid  points:  observed  fluores¬ 
cence  decay  of  1,2-benzanthracene  in  1,2-propanediol  quenched  by 
0.29  M  CBr4.  At  30°  C.  Dashed  line:  attenpt  to  fit  with  G(t)  « 
exp{t/x}  to  minimize  the  sum  of  the  squares  of  the  qeighted  residuals. 
Solid  line:  fit  with  G(t)  =  exp(-at  -  2b/t)  by  iterative  convolution. 
Best  a  and  b  give  R'  ■  7.5  X  and  D.n  =  0.79  x  10“  6  cm2  sec"l. 

Ay 

Figure  5.  Steady-state  experiments  of  fluorescence  quenching  of  1,2- 
benzanthracene  by  CBr4  in  1,2 -propanediol  at  different  temperatures: 

(A)  15°  C;  (B)  25°  C;  (C)  35°  C.  Solid  lines:  calculated  from  Eq. 

(10)  with  parameters  from  Table  1  (our  model).  Dashed  line:  calcu¬ 
lated  from  Eq.  (20)  for  15°  C. 

Figure  6.  (1)  Flash  lamp,  (2)  observed  fluorescence  decay  of  pyrene 

in  mixing  75%  cyclohexanol,  25%  cyclohexane  as  solvent  at  22°  C, 
quenched  by  biacetyl  (0.0677  mole  Jt“l).  (a)  Attempt  to  fit  with 
G(t)°°exp(-t/x)  to  minimize  the  sum  of  the  squares  of  the  weighted 
residuals,  (b)  Fit  with  G(t)  =  exp(-at  -  bt5*)  by  iterative  convolu¬ 
tion  (a  =  6.3  x  107  M"l,  b  =  7.1  x  10^  ns-*5) .  £  channel  =  0.44  ns. 

Figure  7.  (a)  Deconvoluted  curve  of  fluorescence  decay  of  pyrene,  in 

mixing  75%  cyclohexanol,  25%  cyclohexane  as  solvent,  at  22°  C,  quenched 
by  biacetyl  (0.0677  mol  £“*)  by  using  Fast  Fourier  transform  programme, 
(b)  (1)  Flash  lamp;  (2)  observed  fluorescence  decay  of  pyrene; 

(3)  flash  lamp  convoluted  with  the  curve  a.  £  channel  =  0.285  ns. 

Figure  8.  Variations  of  (a)  and  (b)  with  (biacetyl). 

Figure  9.  Steady-state  fluorescence  experiments:  variation  of  <j>0/<j> 
versus  (CBr4) .  (1)  r\  =  1.32  cp;  (2)  n  =  2.04  cp;  (3)  n  =  2.47  cp; 

(4)  n  =  3.02  cp;  (5)  n  =  6.72  cp;  (6)  n  -  11.63  cp;  (7)  n  =  24  cp, 

G  =  22°  C. 
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SECTION  II 

EXCIPLEX  PHOTOPHYSICS 


While  exciplexes  have  been  known  for  several  decades,  there 
has  been  remarkably  little  work  reported  in  which  the  details  of 
the  kinetics  are  examined  with  sufficient  precision  and  over  a 
sufficient  range  of  condition  such  that  a  critical  test  of  proposed 
kinetic  models  is  possible.  Even  the  simple  "Excimer"  model  for 
the  kinetics  involves  a  minimum  of  six  rate  constants.  In  what 
follows  the  recent  work  along  these  lines  done  to  elucidate  the 
details  of  exciplex  photokinetics  is  reports. 


The  photophysics  of  systems  exhibiting  molecular  association 
both  in  the  ground  and  excited  states  has  been  studied.  The  emphasis 
has  been  on  kinetic  models,  the  measurement  of  rate  constants 
associated  with  these  models,  and  the  determination  of  activation 
parameters  and  equilibrium  thermodynamic  parameters  associated  with 
the  exciplex  formation  and  disappearance.  Studies  of  solvent  effects 
and  steric  effects  on  the  behaviour  of  exciplex  systems  have  been 
carried  out.  The  case  of  rapid  equilibrium  where  the  monomer  and 
exciplex  decay  with  the  same  rate  constant  has  also  been  examined. 


1.  INTRODUCTION 


The  terms  "excimer"  and  "exciplex"  are  used  *o  describe  molecu¬ 
lar  or  atomic  aggregates  which  are  unstable  in  their  ground  states 
but  possess  significant  binding  energy  in  their  electronical ly  exci¬ 
ted  states.  The  term  excimer  was  first  used  by  Stevens  end  Hutton1 
to  describe  aggregates  of  like  molecules.  The  term  exciplex  is  due 
to  Walker,  Bender  and  Lurary^  and  pertains  to  the  association  of  unlike 
species.  Mataga,  Okada  and  Ochari^  introduced  the  descriptive  term 
"heteroexcimer"  for  the  exciplex  and  both  are  in  common  use  today. 
However,  the  phenomena  described  by  these  terras  predates  the  nomen¬ 
clature.  Forster  and  Kasper4  first  observed  the  characteristic  red- 
shifted  exciplex  emission  in  concentrated  pyrene  solutions  and  Leon- 
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hardt  and  Weller®  were  the  first  to  observe  exciplex  emission  in  the 
system  perylene  quenched  by  N,N-dimethylaniline. 

Excited  states  of  EDA  (electron  donor- accept  or)  conqplexes  formed 
by  direct  excitation  into  the  CT  band  have  many  properties  in  common 
with  exciplexes  and  in  general  differ  only  in  the  ground  state  sta¬ 
bility.  6»/ 

Excimers  and  exciplexes  play  an  important  role  in  photochemical 
transformations.®  In  the  exciplex,  one  has  exchange  and  charge  trans¬ 
fer  interactions  which  hold  the  pair  in  close  proximity  allowing  for 
the  possibility  of  dimerization.  In  the  exciplex,  one  has  not  only  a 
large  charge  separation  in  the  bound  excited  complex,  but  radical 
ions  can  be  generated  in  the  exciplex  dissociation  process  in  polar 
media.  In  addition,  triplet  formation  via  the  exciplex  is  well  estab¬ 
lished.  Thus,  in  polar  solvents,  one  frequently  sees  radical  ion 
chemistry  from  exciplex  forming  systems  with  the  "photochemistiy" 
acting  merely  to  generate  radicals.  Excited  charge- transfer  complexes 
can  also  serve  as  intermediates  in  photochemical  reactions.9 


2.  PRIMARY  PROCESSES  IN  EXCIPLEX  SYSTEMS 


Two  aspects  related  to  the  primary  processes  in  exciplex  systems 
have  generated  considerable  controversy  over  the  past  15  years.  One 
concerns  the  question  of  radical  ion  generation  directly  from  the 
encounter  complex  prior  to  relaxation  to  form  the  bound  exciplex. 

The  other  question  concerns  the  formation  of  triplet  states  through 
the  donor- acceptor  interaction,  but  without  the  bound  exciplex  acting 
as  an  intermediate  (the  so-called  "fast"  triplet  production) . 

The  first  problem  had  its  origin  in  the  observation10'*1  that 
the  exciplex  emission  yield  fell  off  much  more  rapidly  than  the  exci¬ 
plex  lifetime  as  the  solvent  dielectric  constant  was  raised  from  2.5 
to  10.  Mataga  and  co-workers  initially  attributed  this  to  a  decline 
in  the  exciplex  radiative  lifetime  as  well  as  to  enhanced  radiation¬ 
less  processes  as  the  solvent  polarity  increased. 12 

13 

Weller  and  co-workers  argued  that  as  the  polarity  of  the  sol¬ 
vent  increased,  the  direct  generation  of  radical  ions  took  place  from 
the  encounter  complex.  These  ideas  can  be  combined  in  the  following 
scheme: 


*A*  +  D  (*A* . .  .D) 


(AD)*  -*■  D  +  A 


hvr 


^  C2As"...2Ds+)  -  2 A "  +  V 


2  -  2  + 

where  (  Ag  . . .  D  )  is  the  solvent  shared  ion  pair.  Laser  photocon- 


S.- 
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ductivity  studies  of  Mataga  and  co-workers  provide  evidence  for  the 
direct  generation  of  the  solvent  shared  ion  pair  from  the  encounter 
complex.  However,  they  also  observed  the  growth  of  conductivity 
matching  the  decay  of  the  exciplex  thus  confirming  the  path  whereby 
the  exciplex  thermally  dissociates  into  radical  ions.  There  is  also 
direct  evidence  for  the  sensitivity  of  the  rate  of  exciplex  emission 
to  solvent  polarity,15  but  the  effects  are  not  large.  Thus,  the 
scheme  presented  above  is  probably  correct  in  general,  but  with 
solvent  effects  present  in  most  if  not  all  of  the  individual  steps. 

In  non-polar  solvents,  the  energy  of  dissociation  of  the  exciplex 
is  high  and  the  exciplex  does  not  dissociate  into  radical  ions  to  a 
significant  extent.  This  simplifies  the  photophysics  somewhat. 

The  second  controversial  topic  relates  to  triplet  production 
and  can  be  illustrated  by  two  quite  different  schemes: 

1  l  Tsr  t 

(a)  A*  +  D  +  (AD)*  °A*  +  D 

(b)  *A*  +  C  -*■  [2A"(+)  +  V(+)]  -  [2A'(0  +  2D(+)]  -*■  3A*  +  D 

In  non-polar  solvents,  it  appears  that  process  (a)  is  the  only 
important  one  and  that  process  (b) ,  which  would  generate  the  triplet 
via  a  fast  process  which  bypasses  the  exciplex,  does  not  occur  to  a 
significant  extent.16  However,  in  highly  polar  solvents,  the  weight 
of  the  evidence  points  toward  process  (b)  and  the  existence  of  a 
"fast"  triplet  generation  route.  The  currently  popular  model  involves 
hyperfine  coupling  interactions  in  the  radical  ion  pair  in  highly 
polar  solvents.16,17 


3.  KINETIC  ASPECTS  OF  EXCIJ^R  AND  EXCIPLEX  PHOTOPHYSICS 
The  simples  reaction  scheme  for  quenching  is 

A  -~-V-  -*  A*  +  Q  -*■  (AQ)  *  -*■  products 

3 

A  +  hv  A  or  A 

In  this  scheme,  there  is  no  regeneration  of  A*  from  (AQ)*  and  the 
photokinetics  correspond  to  ordinary  fluorescence  quenching  kinetics. 
(AQ)*,  being  non-emi&sive  in  this  scheme,  remains  an  elusive  tran¬ 
sient  species,  detectable  perhaps  by  absorption  spectroscopy.  Its 
presence  as  an  intermediate  may  be  inferred  from  the  products 
observed  (e.g.,  products  of  radical  ion  reactions),  but  frequently, 
one  can  only  guess  as  to  the  structure  and  properties  of  this  inter¬ 
mediate.  The  rate  of  the  forward  or  quenching  reaction  will  probably 
be  diffusion  controlled  which  limits  the  information  available  from 


V 


40 


photophysical  studies.  Nevertheless,  the  transient  may  be  a  genuine 
exciplex  which  is  merely  very  unstable  relative  to  products  in  the 
solvent  used. 

The  next  level  of  complexity  introduces  reversibility  and  exci¬ 
plex  fluorescence  (hv^) . 


A  +  hv_  A  or  A  A  +  Q  +  hv„ 

F  E 


Whether  or  not  exciplex  emission  is  measurable  depends  on  equipment 
sensitivity,  the  ratio  kg/(k4  +  k5  +  k$) ,  and  the  ratio 

k3[Q]/(k1  +  k2  ♦  k3). 


These  photokinetics  require  that 


-A.t  -A,t 


A.  ~ 

[A*](t)  =  a^  +  a2e 

(1) 

-A,t  -A„t 

C(AQ) *](t)  =  a3(e  +  e  ) 

(2) 

where 

X1>2  =  1/2  [kx  ♦  k2  ♦  k3[Q]  ♦  k4  +  kp  +  {(kj  ♦  k2  ♦ 

♦  k3[Q]  -  *4  -  kp)2  ♦  4k3k4[Q]}5s] 

and 


(3) 

(4) 


Steady-state  analysis  gives  (where  fyq  and  represent  the 
quantum  yields  of  monomer  and  exciplex,  respectively) 


r$M°  k3kp[Q] 

t  7  (k4  *  kp^kl  +  k2-) 


(5) 
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♦g  (k5|  k3[Q] 

TlM  757"rV  (6) 

In  the  above  equations,  it  has  been  assumed  that  kj  j*  f (t) .  Rate 
constants  can  then  be  determined  as  follows: 

k^:  from  the  unquenched  lifetime  (t)q_q  =  (kj  ♦  k2)~*  and  the 
quantum  yield  4ty°  ■  kj/(kj  +  k2)  . 
k2=  from  k^  and  (t)q=0. 

k_:  from  the  slope  of  (X^  +  X2)  vs.  [q],  since  Xj  +  X2  =  (k  + 

k2  +  k3tQ3  +  k4  +  kp*  ,  .  _ 

k.  and  kp:  from  a  plot  of  X]/2  vs  LQJ  since  X^2  =  (kj  + 

4  k2)tk4  ♦  kp)  ♦  kpksCQ].  One  obtains  k4  +  kp  from  the  inter¬ 
cept  and  kp  from  the  slope.  Since  k3  is  known,  k4  is  cal¬ 
culated  from  k4  +  kp.  Problems  obviously  arise  here  if 
k4  «  kp. 

k,.:  from  Eq.  (6),  k5  is  obtained  by  plotting  vs  CqD - 

The  slope  gives 


k,  k.  +  k 
14  p 

from  which  ks  is  calculated  since  all  other  rate  constants 
are  known. 

kt :  since  k  =  kc  +  k,,  once  kr  is  known,  one  can  calculate  k, . 
o  p  b  o  5  o 

Alternate  methods  are  discussed  elsewhere  for  evaluating  some 
of  the  rate  constants  in  the  above  exciplex  photokinetic  scheme,  but 
they  appear  to  offer  no  real  advantage  over  the  direct  approach  out¬ 
lined  above. 

Once  the  constants  (kj  +  k2)»  k-j,  k4  and  kp  have  been  determined, 
one  can  calculate  the  Stem-Volmer  constant  Kgy 


IqJ  =  (kj  +  k2)(k4  ♦  k  ) 


and  compare  with  experimental  values  derived  from  steady-state  mea¬ 
surements.  Discrepancies  are  interpreted^»20  in  terms  of  failure  of 
the  assumption  that  kj  f  f(t)  or  a  failure  of  the  kinetic  scheme  to 
describe  the  system  in  question. 

Rate  constants  may  then  be  determined  as  a  function  of  tempera¬ 
ture.  When  linear  In  k^  vs  1/T  plots  are  obtained,  one  can  calculate 


'J! 
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the  parameters  in  the  equation 

-ae//rt 

k.  =  A.e  (8) 

From  the  temperature  coefficient  of  the  ratio  of  one  obtains 

AHC°,  the  enthalpy  of  the  excited  state  equilibrium  and  from  A3  and 
A4  one  may  obtain  the  entropies  of  activation  of  the  forward  and 
reverse  reactions  in  the  excited  state  equilibrium.  Finally,  from 
energy  conservation  (see  Fig.  1) 

hv(A  -  A*)  =  AHC°  +  (hVg)^  +  AHr  (9) 


where  (hVE)max  is  the  maximum  in  the  exciplex  emission  and  AHr  is  the 
repulsive  potential  energy  in  the  Franck-Condon  ground  state  at  the 
equilibrium  separation  of  A  and  Q  in  the  exciplex.  This  yields  AHR. 

It  is  emphasized  that  the  rate  constants  IC3,  k4  and  kp  are 
obtained  without  recourse  to  steady-state  data.  It  is  only  when  one 
calculates  ki,  k2,  kg  and  kg  that  it  is  necessary  to  use  4^°  and 
4>e/4>u*  the  latter  ideally  as  a  function  of  [Q]. 

Likewise,  one  obtains  AHC°  without  recourse  to  the  arguments  of 
Stevens  and  Ban,21  i.e.,  if  k4  »  kp,  then 


(10) 


If  this  is  true,  then  of  course  one  obtains  AHC°  from  In  (^/(^[Q])  ^ 
vs  1/T  at  constant  [Q] .  As  Selinger  and  McDonald22  point  out,  there 
are  many  cases  where  it  is  not  clear  that  this  is  a  valid  approxima¬ 
tion.  One  obviously  needs  individual  rate  parameters  to  establish 
AHC°  with  certainty,  unless  one  can  measure  4’e/4>M  °ver  so  wide  a 
tenperature  range  that  one  sees  the  limiting  behaviour  for  k4  »  kp 
as  well  as  for  k4  «  kp.  Finally,  it  is  emphasized  that  Kgy  is  a 
function  of  k3,  k4,  ks  and  kg  as  well  as  (t)q=q.  Changes  in  Kgy  are 
clearly  almost  inpossible  to  inteipret  unless  one  has  values  for  the 
individual  rate  constants,  since  all  are  potentially  functions  of 
both  temperature  and  solvent!  Kgy  is  in  fact  not  a  particularly 
informative  parameter  except  as  a  route  to  individual  rate  constants. 


Some  of  the  reasons  for  believing  that  the  above  model  is  real¬ 
istic  for  many  exciplex  systems  in  non-polar  or  slightly  polar  sol¬ 
vents  are  as  follows:  fa)  Two  conponent  fluorescence  decay  predicted 
’  y  Eq.  (1)  is  observed, *8, 19  along  with  the  growth  and  decay 
predicted  in  Eq.  (2).  (b)  Rate  constants  derived  from  transient 

luminescent  measurement  agree  with  steady-state  measurements  to  the 
extent  expected  from  the  theory  of  diffusion  controlled  reactions.1® 


'mi 


(c)  The  value  for  kj  is  reasonable  -  generally  near  but  never  above 
the  diffusion  control  limit,  (d)  The  observed  variation  of  Xj.and 
X2  with  [Q]  can  be  reproduced  by  a  single  set  of  rate  constants. 

Thus,  while  one  of  course  cannot  rule  out  more  conplex  schemes, 
we  have  found  no  need  to  envoke  them  to  explain  numerous  experimental 
observations  in  non-polar  or  slightly  polar  solvents. 

Extensive  studies  over  the  last  15  years  by  a  number  of  groups 
indicate  that  the  same  general  conclusion  is  valid  in  excimer  systems. 


The  process 

(AQ)  *  +  Q  -*■  Products 

23 

appears  important  only  at  high  concentrations  of  quencher,  but 
most  experimental  studies  have  been  below  the  quencher  level  where 
this  is  kinetically  significant. 

The  complication  represented  by  the  process 

A  +  Q  v .  ^  (AQ)  (AQ)  * 

is  more  difficult  to  rule  out.  While  the  exciplex  is  by  definition 
repulsive  in  the  ground  state,  some  systems  are  certain  to  represent 
borderline  cases  of  weak  attraction.  Such  weak  attraction  causes  the 
distribution  of  Q  around  A  to  be  non-random.  This  alters  the  kinetics 
of  the  forward  reaction  (k 3) ,  which  are  already  complex  due  to  the 
time  evolution  of  concentration  gradients jn  addition,  one 
must  consider  the  photoselection  of  AQ  pairs  if  one  is  using  an 
exciting  wavelength  which  favours  AQ  over  A. 

Unfortunately,  one  does  not  always  see  two  components  in  the 
fluorescence  decay  of  a  fluorophor,  even  when  the  emission  spectrum 
clearly  indicates  emission  from  both  A*  and  (AQ)*.  Measurements  in 
the  spectral  region  in  which  the  decay  is  expected  to  be  predominately 
due  to  the  exciplex  are  frequently  found  to  give  similar  if  not  iden¬ 
tical  decay  curves  as  obtained  from  measurements  of  the  decay  of 
quenched  A*.  In  many  cases,  the  decay  is  very  close  to  a  single 
exponential  with  a  single  decay  constant  X,  and  frequently  X  is 
found  to  be  a  strong  function  of  [Q].  It  is  not  unusual  to  find  for 
such  systems  that  X  follows  the  equation 

X  =  (Xq  +  B[Q])/(1  +  C[Q])  (11) 

If  the  following  approximation  is  valid 
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C(k4  ♦  k3£Q])2  ♦  2(k4  -  k3lQ])Ckp  -  kj  -  kjY 
2(k.  -  k,[Q])Ck  -  k.  -  k  ) 

~  i.  .1.  .43  p  1  2 


=  k4  +  ks[Q] 


k,  +  kJQJ 


Eq.  (3)  can  be  reduced  to  Eq.  (11)-  This  is  related  to  the  well 
known  condition 

k4,  k3[Q]  »  kj,  k2,  kp  (] 

to  which  must  be  added  the  more  restrictive  condition 

k4  +  k3[Q]  >  2(k4  -  k3[Q])Ckp  -  kj  -  k2)  (1 

The  conditions  quoted  in  Eqs.  12  and  14  thus  yield 

,  »„♦ 


- 1  *  Ck3/k4)LQj 

A  plot  of  (A  -  A0)  1  vs  1/CQ]  should  yield  a  straight  line  with  slope 
S  given  and  intercept  I  given  by 

S'1  ■  Ckp  -  i0) (k3/k4) 

(16) 

I"1  =  k_  -  A 


This,  S/I  =  k4/k?  and  I  +  =  The  [Q]  dependence  alone  of  A 

(along  with  AQ)  fails  to  yield  the  complete  set  of  rate  constants, 
but  rather  yields  only  kp  and  ratio  k3/k4.  However,  from  the  steady- 
state  quenching  measurements,  we  have 


i/k3  ■  i/(kq)ss 


where 


<V«s  -  VT0  (18) 

When  this  equation  is  used  with  values  of  k3/k4  and  kp,  one  can  obtain 
kj.  Then  k4  can  be  deduced.  It  is  important  to  recognize  that  the 
success  of  this  approach  is  the  analysis  of  the  so-called  rapid 
equilibrium  case  depends  critically  upon  being  able  to  obtain  a 
significant  number  from  the  subtraction  of  (l/kq)ss  and  k4/k3kp.  It 
is  not  uncommon  that  the  two  terms  in  Eq.  (17)  are  nearly  identical. 
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which  causes  the  analysis  to  fail. 

Figures  2,  3  and  4  illustrate  experimental  data  consistent  with 
Eqs.  (15).  and. (16).  The  use  of  steady-state  (Fig.  5)  data  with 
Eqs.  (17)  and  (18)  yielded  individual  rate  constants.1®  As  will  be 
discussed  below,  by  varying  the  temperature,  one  can  cause  the  cyano- 
naphthalene-olefin  system  to  exhibit  the  classical  two  component 
decay.  The  rate  constants  obtained  are  consistent  with  those  obtained 
in  the  rapid  equilibrium  region. 18 


3.  EFFECT  OF  SOLVENT  AND  TEWERATURE  ON  PH0T0K1NETICS 

Consider  first  a-cyanonaphthalene  (CCN)  quenched  by  1, 2-dimethyl 
cyclopentene  (DMCP) .  When  this  system  is  examined  in  the  vapour 
phase, 24  the  quenching  is  observed  to  be  very  inefficient.  One  can 
estimate  from  lifetime  measurements  that  the  effective  rate  constant 
for  quenching  is  about  10" ^  of  the  gas  kinetic  limiting  rate.  One 
sees  single  exponential  decay  provided  a  buffer  gas  is  present,  but 
the  system  has  thus  far  resisted  attempts  to  recover  rate  constants. 

1 8 

If  one  now  examines  this  same  system  in  hexane,  one  finds  that 
the  observed  photokinetic  behaviour  is  a  strong  function  of  tempera¬ 
ture.  At  room  temperature  and  above,  this  system  provides  a  classic 
example  of  the  rapid  equilibrium  case.  Eq.  (15)  is  followed  exact'y 
and  rate  constants  can  be  derived. 18  The  fluorescence  decay  follows 
a  single  exponential.  If  one  now  decreases  the  temperature  at  about 
0°  C,  two  component  fluorescent  decay  becomes  measurable  (see  Figs. 

6  and  7)  and  rate  constants  can  be  obtained  by  the  analysis  described 
above.  The  agreement  between  calculated  and  observed  values  of  Xj 
and  X2  is  shown  in  Figs.  8  and  9.  These  rate  constants  are  consistent 
with  those  obtained  from  the  analysis  of  the  rapid  equilibrium  case.18 
The  temperature  range  of  40°  to  -40°  C  takes  one  from  one  limiting 
case  to  the  other.  The  reason  for  this  is  that  k4  varies  by  three 
orders  of  magnitude  over  this  tenperature  range.  At  a  low  tempera¬ 
ture,  the  reverse  reaction  is  slowed  down  and  one  can  see  two  compo¬ 
nent  decay,  whereas  at  high  temperature,  the  excited  state  quasi¬ 
equilibrium  is  the  dominant  kinetic  process  and  the  other  processes 
merely  bleed  off  excited  species  from  this  quasi-equilibrium.  This 
is  illustrated  in  Fig.  10. 

Not  only  does  k4  control  the  photokinetic  behaviour  in  non-polar 
solvents,  but  also,  it  rather  than  kp,  is  the  most  solvent-sensitive 
rate  constant  in  this  system. 25  Merely  changing  the  solvent  from 
hexane  to  diethyl  ether  decreases  k4  at  a  given  tenperature  by  a 
factor  of  20.  If  one  uses  a  solvent  such  as  THF,  the  photokinetics 
are  now  simply  those  of  a  diffusion  controlled  reaction.  The  solvent 
polarity  has  presumably  reduced  k4  to  the  point  where  there  is  no  sig¬ 
nificant  reversibility,  and  individual  rate  constants  cannot  be 
recovered  by  nanosecond  fluorometry. 
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A  study  ’  of  the  tenperature  coefficients  of  IC3  and  IC4 
revealed  that  AHC°,  the  enthalpy  change  of  the  excited  state  reac¬ 
tion,  is  independent  of  the  solvent  change  from  hexane  to  ethyl 
acetate,  as  is  AE^,  the  activation  energy  for  the  reverse  reaction. 
The  dramatic  drop  in  k4  with  solvent  is  thus  due  to  the  pre-exponen¬ 
tial  term  and  can  be  interpreted  as  an  entropy  of  activation  effect 
associated  with  changes  in  specific  solvation  by  these  two  solvents 
of  the  highly  dipolar  exciplex. 

IC4  also  appears  to  depend  upon  the  ionization  potential  of  the 
donor  when  COJ  is  quenched  by  various  olefins  and  a  diene  in  hexane. 
As  the  ionization  potential  of  the  donor  decreases,  IC4  dramatically 
increases.1 8  This  phenomenon  deserves  much  more  detailed  study 
including  the  recovery  of  thermodynamic  parameters  for  a  number  of 
donors . 


kp  is  also  observed  to  decrease  as  one  increases  the  solvent 
polarity.  This  is  in  part  due  to  a  drop  in  kg,  but  the  largest 
contribution  is  a  drop  in  kg. 


26-  2  8 

Two  products  are  observed  by  other  laboratories  ”  for  the 
reaction  of  excited  CCN  with  DMCP.  These  are  a  ring  adduct  (I)  and 
a  nitrile  adduct  (II).  Lewis  and  Hoyle^?  have  in  fact  observed  that 
nitrile  addition  is  favoured  in  hexane  and  ring  addition  in  ethyl 
acetate  and  that  the  quantum  yield  of  each  product  decreases  with 
increasing  solvent  polarity.  This  is  consistent  with  our  results 
for  the  variation  of  kg  with  solvent  polarity  if  one  assumes  that  k6 
is  predominantly  the  product  formation  rate  constant. 


A  similar  trend  is  found  in  the  kg  values,  kg,  which  was  cal¬ 
culated  with  Eq.  (5)  on  the  assumption  that  kj,  the  CNN  fluorescence 
probability,  is  tenperature  independent,  decreases  with  increased 
solvent  polarity.  Owing  to  the  problem  of  evaluating  the  contribu¬ 
tions  of  various  component  configurations  to  the  total  exciplex  wave 
function,  a  rigorous  calculation  of  the  radiative  transition  probabil¬ 
ity  is  not  possible.  However,  the  observed  decrease  in  kg  supports 
the  opinion  that  the  electronic  structure  of  the  complex  changes  as 
the  solvent  polarity  is  changed,  k  is  related  to  the  transition 


(I) 


(II) 


47 


dipole  moment,  for  the  exciplex  through  the  equation 


29 


*E 


3hc' 


,,  4  3 

64  it  nv 

e 


(19) 


Mg  in  turn  is  given  by 

% 


^E^op^G* 


(20) 


in  which  Mop  is  the  electric  dipole  moment  operator,  and  ip  and 
are  given  by 

ip£  =  a^fAV)  ♦  a24<2(A+D‘)  +  a3<|>3(A*D)  +  a^(AD*)  (21) 

iPG  =  iPo(AD)  (22) 

Because  the  ground  state  is  dissociative,  it  has  been  assumed  for  the 
present  argument  that  charge-transfer  contributions  to  its  wave  func¬ 
tion  are^negligible.  If  the  second  and  fourth  terms  in  are  neg¬ 
lected,  M  is  given  by 

Mg  =  ai  <  iKaV)|m  |*o(AD)  >  -a3  <  *3(A*D)  |y 0  [i|»0(AD)>  (23) 


OP 
-  + 


The  matrix  element  <^(A  D  ) | M0p I 'J'o (AD) >  is  approximately  proportional 
to  the  overlap  of  the  donor  orbital  in  D  and  the  acceptor  orbital  in 
A; 29  consequently,  it  is  rather  small  and  decreases  with  increasing 
intermolecular  distance.  As  the  exciple^c  becomes  more  ^onic,  a  3 
decreases;  thus,  the  dominating  term  in  a$<ipJ[A*D)  |  Mop  I 'J’o  (AD)  >, 

decreases  in  more  polar  solvents.  According  to  this  analysis,  there¬ 
fore,  both  the  increasing  charge -transfer  contribution,  which  reduces 
the  second  term  in  Eq.  (23)  and  the  probable  increase  in  the  inter¬ 
molecular  distance  which  reduces  the  first  term,  contribute  to  the 
decrease  in  k,.  with  increasing  solvent  polarity .  30, 31 

Although  kg  decreases  in  the  more  polar  solvents,  the  quantum 
yield  of  exciplex  emission,  4>_,  given  by 

E 

h '  k„<k4  *  V  *  VjW 

or,  at  infinite  quencher  concentration,  by 
00  ^5 

*e  -  r 


(24) 


(25) 


increases.  In  general,  in  exciplexes,  decreases  with  increasing 
solvent  polarity; 32  an  increased  has  also  been  observed  in 
media  of  moderate  dielectric  constant. 33  Obviously,  when  increased 
solvent  polarity  induces  an  increase  in  kp  or  a  high  probability  for 
instantaneous  radical  ion  production,  with  rate  constant  k.3’,  4>£  will 
decrease.  4>_  is  then  given  by 


k5k3 


(k,  .  k.'Jk  *  k.’k 


3  4 


(26) 


However,  the  present  results  show  that  neither  of  these  effects  occurs 
for  the  CNN-DMCP  system  in  the  studied  solvents.  It  is  worthy  of  note 
that,  in  the  present  range  of  solvent  polarity,  both  increases30*34-  b 
and  decreases3!* 33  in  have  been  observed. 

P 
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We  have  recently  examined  the  quenching  of  pyrene  with  N,N- 
dimethylaniline  (DMA)  substituted  on  the  ring  with  tert -butyl  groups. 
This  work  was  suggested  by  the  observation  of  Taylor  et  al.3®  that 
the  Stem-Volmer  constant  for  quenching  strongly  increases  upon  3,5 
di-tert  butyl  (DBDMA)  substitution.  Hie  question  then  was:  which 
rate  constants  in  the  expression 


k,k 

k  =— Lr 


q  k4  +  kr 


(27) 


were  being  influenced? 

Detailed  kinetic  studies  using  the  method  of  analysis  described 
above  conclusively  demonstrate  that  the  large  increase  in  kq  for  this 
system  in  non-polar  solvents  is  due  to  a  large  decrease  in  k4.  The 
forward  reaction  is  diffusion  controlled  for  both  DMA  and  DBDMA 
quenchers.  kp  decreased  somewhat  upon  di-tert  butyl  substitution. 
These  results  point  to  the  danger  associated  associated  with  inter¬ 
pretation  of  kq.  In  this  case,  one  has  several  remarkable  effects 
occurring  simultaneously.  Di-tert  butyl  substitution  increases  the 
stability  of  the  exciplex  while  altering  the  close  sandwich  pair 
configuration.  Thus,  there  are  both  electronic  and  steric  effects  in 
operation  and  considerable  information  can  be  obtained  from  detailed 
kinetic  studies  as  a  function  of  solvent  and  temperature.37  kq  alone 
tells  one  almost  nothing! 


The  formation  of  exciplexes  and  excimers  is  almost  always  at  the 
diffusion  controlled  limit.  Thus,  the  formalism  presented  above,  and 
almost  universally  used  to  analyze  the  photokinetics  of  these  systems, 
is  an  approximation  to  the  truth  because  in  fact  k3  is  a  complex  func¬ 
tion  of  time.  This  subject  has  received  detailed  attention  in  another 
lecture  by  the  author,  but  the  subject  requires  some  conment  in  the 
present  context. 
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Mention  must  first  be  made  of  the  studies  of  Eisenthal  and 
co-workers."*9  They  employed  a  mode-locked  laser  for  picosecond  exci¬ 
tation  and  were  able  by  optical  absorption  to  observe  the  growth  of 
the  exciplex  between  anthracene  and  diethylaniline  in  hexane.  Their 
data  fit  the  theory  of  diffusion  controlled  reactions  quite  well  at 
the  level  of  sophistication  described  by  Smolochowski40  including 
terms  normally  ignored,  but  important  at  short  times. 

19 

Ware  and  Hui  have  also  examined  this  problem  by  introducing  an 
approximate  time  dependent  for  kj  into  the  differential  equations 
for  the  exciplex-excimer  kinetic  scheme,  i.e., 

4j£l  -  k4[*Q*]  -  [kj  *  k2  .  k3(t)[Q]][»*] 

«  kjCt)[Q][A.]  -  [k4  ♦  k5  .  k6][AQ*] 

This  system  of  equations  was  then  solved  by  numerical  integration. 

It  was  found  that  distortions  to  the  leading  edge  of  both  the  exciplex 
growth  curve  and  the  monomer  decay  curve  were  caused  by  the  time 
dependence  of  k^. 

However,  if  deconvolution  was  carried  out  only  over  those  por¬ 
tions  of  the  decay  curves  where  kj(t)  had  reached  its  limiting  time 
independent  value,  the  problem  could  be  avoided. *9 


4.  FLUORESCENCE  QUENQIING  VIA  H-BOND  FORMATION 


The  techniques  described  above  are  readily  applicable  to  the 
quenching  of  excited  species  by  H-bond  formation.  However,  one 
frequently  has  the  added  complication  of  a  bound  ground  state  complex 
which  competes  with  A*  for  the  incident  radiation.  The  photokinetic 
scheme  can  now  be  depicted  by 


AH  +  Q 


I 


abs 


x<5 


AH*  +  Q 


To  -  (k1+k2) 


AH— Q 


Ijabs*1^ 
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A  number  of  years  ago,  the  author  piilished  an  account  of  the  behav¬ 
iour  of  such  a  system  where  kj  >>  k^.  Recently,  Dr.  M.  Martin  and 
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and  the  author  have  completed  a  study  in  this  laboratory  of  a  more 
difficult  system,  carbazole  quenched  by  pyridine  and  siistituted 
pyridines.  This  system  exhibits  weak  fluorescence  from  the  excited 
H-bond  complex.  Steady-state  analysis  of  the  above  scheme  yields 


kT'(l  -  <5)  k5T' Cl  -  5) 

5  +  — ; - ...  —  + 


■  t. 


1  ♦  k-T1 
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x  k,  \ 


^V’.0-45) 

V  k3 

\  1  *  v  I 

ksV'd  - «) 


k,(l  •  k,!’)2' 


[Q] 


(28) 


whereas  the  transient  behaviour  is  given  by  modifications  of  Eqs. 

(1)  and  (2),  provided  the  excitation  is  almost  entirely  to  A*  for 
the  lifetime  measurements.  In  this  system,  the  absorption  due  to  the 
ground  state  complex  is  somewhat  separated  from  that  of  the  carbazole. 
Our  kinetic  studies  have  resulted  in  the  following  conclusions: 

(a)  The  kinetic  model  given  above  appears  valid.  Rate  cons¬ 
tants  derived  from  the  kinetics  by  transient  measurements,  when  intro¬ 
duced  into  Eq.  (28),  yield  excellent  fits  such  as  that  illustrated 

in  Figs.  11  and  12.  The  only  adjustable  parameter  is  the  excited 
state  complex  lifetime  x'  which  is  found  by  an  iterative  procedure. 

For  the  carbazol e-pyridine  excited  complex,  we  obtain  a  lifetime  of 
52  psec.  Thus,  there  are  very  rapid  radiationless  processes  in  this 
excited  state  complex  which  probably  yield  the  ground  state  molec¬ 
ule.  ^  No  photochemistry  or  proton  transfer  is  observed. 

(b)  The  effect  of  methyl  and  phenyl  s;i>stitution  on  the  pyridine 
has  also  been  studied  and  leads  to  the  conclusion  that  the  coplanarity 
of  the  carbazole  and  pyridine  ring  systems  is  not  significant  in 
determining  the  quenching  probability. 

(c)  The  deuteration  of  the  amine  hydrogen  of  the  carbazole  also 
produces  no  change  in  the  photophysical  behaviour,  which  suggests  that 
N-H  Franck-Condon  overlap  is  not  involved  in  determining  the  radia¬ 
tionless  process  in  this  excited  complex.  Thus,  the  notion^  that 
quenching  involves  a  delocalization  type  interaction  between  the  it 
systems  of  the  carbazole  and  pyridine  is  probably  not  correct. 
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5.  OTHER  TOPICS 

Time  and  space  limitations  obviously  make  it  inpossible  to  deal 
with  all  aspects  of  excimer-exciplex  photophysics.  The  emphasis  in 
in  this  work  has  been  on  the  determination  of  rate  constants  by 
pulse  fluorometry  and  steady-state  fluorescence  spectroscopy  and 
relative  steady-state  intensity  measurements. 

For  more  information,  the  reader  should  consult  the  reviews 
by  Mataga  and  Ottolenghi, Davidson,®  Ware,44  Stevens,4®  Birks,4^ 
Ware  and  Gordon, 4?  Beens  and  Weller,®  as  well  as  the  classical  works 
of  Birks4®  and  Mataga  and  Kubota.4® 

Examination  of  these  reviews  will  reveal  the  great  inpact  of 
laser  spectroscopic  techniques  in  this  field  as  an  approach  to  under¬ 
standing  the  structure  of  the  excimer  and  exciplex  which  compliments 
pulse  fluorometry  studies. 


6.  CONCLUSIONS 

The  molecular  conplex  stable  in  its  ground  electronic  state 
has  been  studied  extensively  over  the  past  several  decades  and  while 
there  are  many  unaswered  questions,  the  properties  and  behaviour  of 
such  species  is  reasonably  well  understood.  The  exciplex  offers  a 
much  greater  challenge.  One  cannot  study  its  physical  and  chemical 
properties  at  one's  leisure.  It  cannot  be  crystallized  and  subjected 
to  x-ray  analysis.  One  must  depend  upon  indirect  observations,  the 
interpretation  of  which  is  rarely  straightforward.  It  is  reasonable 
to  expect  that  as  we  learn  more  about  the  physical  and  chemical 
properties  of  the  exciplex,  we  will  see  many  more  applications  of 
this  particular  electron  donor-acceptor  phenomenon.  One,  however, 
should  not  underestimate  the  time  required  to  collect  quantitative 
kinetic  and  thermodynamic  data.  To  study  a  single  system  as  a  func¬ 
tion  of  solvent,  temperature  and  concentration  of  quencher  requires 
one  half  of  a  man  year,  provided  there  are  no  unexpected  problems  and 
the  system  is  "cooperative".  If  one  so  wishes  to  study  in  this  detail 
a  whole  set  of  donors  and  acceptors,  the  magnitude  of  the  task  is 
somewhat  frightening  as  well  as  expensive.  Add  to  this  studies  of 
photochemical  products  and  laser  flash  photolysis  studies  of  tran¬ 
sient  spectra,  and  one  may  begin  to  visualize  the  magnitude  of  the 
exciplex  problem.  Finally,  all  the  work  should  be  done  in  four  or 
five  solvents. 
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Figure  1.  Schematic  potential  energy  diagram  for  exciplex  forma¬ 
tion.  Eo-0:  0-0  molecular  fluorescence  energy;  AH0:  exciplex 

binding  energy;  h\^:  energy  of  exciplex  fluorescence  maximum; 

Er:  ground  state  repulsion  energy. 

Figure  2.  Typical  exciplex  fluorescence  decay  curves  for  a-cyano- 
naphthalene  unquenched  and  quenched  1, 2-dimethyl cyclopentene  in  hex¬ 
ane  at  25°. 


Figure  3.  Exciplex  fluorescence  decay  constant  in  hexane  as  a  func¬ 
tion  of  1,2-dimethyl cyclopentene  concentration.  Points,  experimental. 
Lines,  best  fit  using  Eq.  (25). 


Figure  4.  Linearization  of  the  curves  shown  in  Fig.  3. 

Figure  5.  Steady-state  Stern- Volmer  plots  for  a-cyanonaphthalene 
quenched  by  1,2- dimethyl cyclopentene  in  hexane. 


Figure  6.  Typical  monomer  decay  curve:  a-cyanonaphthalene  + 
dimethylcyclopentene-1, 2  (2.85  x  10“2  M)  at  -31°  C. 


Figure  7.  Typical  exciplex  decay  curve:  a-cyanoanpthalene  + 
dimethylcyclopentene-1, 2  (1.92  x  10"2  M)  at  -40°  C. 


Figure  8.  Comparison  of  experimental  and  calculated  X^. 

Figure  9.  Comparison  of  experimental  and  calculated  X£. 

Figure  10.  Arrhenius  plots  of  k^,  k^,  k  .  Open  circle  values  were 
taken  from  paper  I.  " 


Figure  11.  Best  fit  obtained  between  experimental  and  calculated 
data  for  the  change  in  the  total  carbazole  fluorescence  yield  (free 
and  bonded  carbazole),  with  increasing  quencher  concentration,  in 
oxygen- free  solution:  (•)  Q  ■  pyridine;  (O )  Q  =  2,6-DMP. 


Figure  12.  Best  fit  obtained  between  experimental  and  calculated 
data  for  the  change  in  the  total  carbazole  fluorescence  yield  (free 
and  bonded  carbazole),  with  increasing  quencher  concentration,  in 
air-saturated  solution:  (£)  Q  ■  pyridine;  d  b  Q  *  2,6-DM’. 
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